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Abstract. With a special 'Ansatz' we analyse the regularity pro- 
perties of atomic electron wavefunctions and electron densities. In 
particular we prove an a priori estimate, swp veB ( xR \ |VV>(y)| < 
C(R) sup yeB ( x 2R) \ip{y)\ and obtain for the spherically averaged 
electron density, p(r), that p"(0) exists and is non-negative. 



Avec un 'Ansatz 7 special nous analysons les proprieties de regularity 
des fonctions d'ondes atomiques et des densites d' electron. En parti- 
culier nous prouvons une estimation a priori, 

sup \ViP(y)\<C(R) sup \iP(y)\ 

y£B(x,R) yeB(x,2R) 

et obtient pour la densite d' electron moyennee sur la sphere p(r), que 
p"(0) existe et est non-negative. 

1. Introduction and Results 

Let V be the Coulomb potential for an atom consisting of a nucleus 
of charge Z (fixed at the origin) and N electrons: 

- Z 1 

V(x) =V(X U ... ,X N ) =J2 ~ ^+ V | r, 

3=1 1 l<j<k<N K| 

x = (xi, ... ,x N ) e R 3N , Xj = (xj^Xj^, x j>3 ) e JR 3 , j = 1, . . . , N, 

(l.l) 

and let H be the corresponding N - electron Hamilton operator: 

H = H N = -A + V (1.2) 

with 

3 d2 



being the kinetic energy operator of the N electrons. The quadratic 
form domain of H is W 1,2 (M. 3N ), see Reed and Simon [Ol. Assume 
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ij; G L 2 (R 3N ) is a real- valued normalised eigenfunction of the operator 
H: 

(H-E)1> = , W\\ = W\\vcr»») = 1- (1-3) 



It is known that then ip is continuous with bounded derivatives, and 
ip G W 2 ' 2 (M. 3N ) (Kato ||) and that ip is in fact analytic away from the 
sine; ularities (in R 3N ) of V, since V is here real analytic (see Hopf M). 
In this paper we derive various qualitative and quantitative properties 
of the wave function ip and of the corresponding one-electron density 



p(x) = / \^{ x i x 2i ■ ■ ■ > xn)\ 2 dx 2 ■ ■ ■ dxN , x G M 3 , (1.4) 

7R3(JV-X) 

as well as of its spherical average (x = ru,r = \x\,u = x/\x\ G § 2 ) 
p(r) = [ p{ru) du 

\if)(ru, x 2 , ■ ■ ■ , xn)\ 2 dx 2 ■ ■ ■ dxjsi dui , r 6 [0,oo). 



§2 J R 3(N-1) 

(1.5) 

Remark 1.1. Aside from Kato's classical results (see Kato the 
local behaviour of electron wavefunctions has been investigated more 
recently by Hoffmann-Ostenhof et al. j|, [§]. The electron density itself 
has been studied extensively in the large- Z-limit, see Lieb 0. Except 
for the spatial asymptotics, see Ahlrichs et al. [|IJ, there are virtually 
no recent rigorous results on p despite the fact that the density is the 
central object in various popular numerical approximation schemes, as 
Density Functional Theory (DFT) and all the various descendants of 
Hartree-Fock theory. 

We now present our results. 

Theorem 1.2. Let ip be as in . For all R G (0,oo), there exists 
a constant C = C(R) such that 

sup |W>(y)| < C sup |V(y)| for all x G M. 3N . 

y£B(x,R) yeB(x,2R) 



Remark 1.3. This result complements the result by Simon [14], Thm. 
C.2.5 (C14)] for the case of operators of the form (|1.2|) , but with V in 
the Kato-class iY n ' 1 (l n ): for 5 G [0, 2) (5 = : n > 3), 

V G K n ' 5 {W l ) <=> lim sup / . ^yj 9 ,. dy = 0. 



J \ x -y\ <e \x-y\ n - 2+d 

The Coulomb potential (1.0) is in K 3N < S (R 3N ) for all 5 G [0, 1), but is 
not in K 3N ^(R 3N ). 

We recall the definition of Holder continuity: 
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Definition 1.4. For C 1" an open set, 6 N, and a G (0, 1], we 
say that the function u belongs to C^"(Q) whenever u G C k (Q), and 
for all (3 G N n with \(3\ = k, and all open balls B(xo,r) C f2, we have 

sup 1 '— , a yUn <C{x ,r). 



As a consequence of the proof of Theorem |1.2| we get: 
Proposition 1.5. Let 

N 

2^ ^Xj-Xk]. (1.6) 

l<j<k<N 





Z 






1 




n in 


(0) 


i> = 


e F cj> 



with 

<\> G Cl^(R 3N ) for all a G (0, 1). 

Remark 1.6. This result classifies the singularities of as those 
coming from V-F: = ipV F + e F V0. Kato || proved that V?/> 

is bounded, but as the ground state of Hydrogen-like systems (N = 
1,E = —Z 2 /4, ijj(x) = c e~ z \ x ^ 2 ,x G M 3 ) shows, it is not in general 
continuous. 



Remark 1.7. The results of Theorem [L2] and Proposition |1.5| easily 
generalise to the case of molecules: L nuclei, of charge Zi, fixed at 
Ri G M 3 , I = 1, . . . , L, with corresponding N - electron Hamilton 
operator 

N / L 

H N ' L = J2 -A, 




We assume throughout when studying p and p that E and r/> in (|L 
are such that there exist constants Co, 7 > such that 

|V>(x)| < C e~ 7|x| for all x G M 3Ar . (1.7) 

For references on the exponential decay of eigenfunctions, see e. g. 
Simon [14]. 

Remark 1.8. Inequality ( |1.7| ) holds when E < inf a ess (H N ). In this 
case, we let e = Eq" 1 — E with E$~ l the ground state energy of the 
(N — 1) - electron operator: 

HN ~ X = E f- A i - A) + E 1 — ~ — r (1-8) 
^ V / ^ Wj - x k \ 

j =2 \ I 3\/ 2<j<k<N 1 3 ft| 
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By the HVZ-theorem (see Cycon et al. [|, Theorem 3.7], inf a ess (H N ) = 
Eq" 1 , and so e > if E < inf a ess (H N ). When we study H in a sym- 
metry sector, so that if> transforms according to this symmetry, then 
Eq" 1 stands for the groundstate energy of the ionized particle system 
described by the Hamiltonian H 1 ^' 1 in the appropriate symmetry sub- 
space as determined by the symmetry behaviour of if). For this case 
a modified version of the HVZ-theorem holds (see e. g. Reed and Si- 
mon [|l^, Thm. XIII. 17'] and Zhislin and Sigalov ]15[]). This includes 
in particular the physically important case of real atoms (Pauli prin- 
ciple). So if E lies below the beginning of the essential spectrum of 
H considered in a symmetry sector, then analogously to the above the 
ionisation energy e > and if) satisfies (|1.7| ). 



Remark 1.9. When assuming ( |1.7j) , Theorem |TT2"| implies that [V^x)) 



also decays exponentially for |x| — > oo. 



Remark 1.10. From (|1.7| ) and Lebesgue's Dominated Convergence 
Theorem follows that the density p is continuous in R 3 . 



Theorem 1.11. Let if> be given according to (L3) and assume that 
(Ol) holds. Then: 



(i) The function p defined in ( |1.4| ) satisfies, in the distributional 
sense, the equation 

1 Z 

--Ap p +h = in R 3 , (1.9) 

2 r 

where 

h e C a {R 3 \ {0}) n L°°(R 3 ) for all a G (0, 1) 
and 

p e C 2 ' a (R 3 \ {0}) n C°'\R 3 ) for all a G (0, 1). 

(ii) The function p defined in ( |1.5|) satisfies 

1 Z ~ 

- -Ap p + h = 0forr>0 (1.10) 

2 r 

where h(r) = j s2 h(ru) duo. Thereby, 

h G C a ((0, oo)) n C°([0, oo)) for all a G (0, 1) 

and 

p G C 2 ' a ((0, oo)) n C 2 ([0, oo)) for all a G (0, 1). 
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(iii) 

h(x) < C{R) ( / p(y)dy + p(x) ) /oradiel 3 , (1.11) 

\Jb(x,R) J 

h(x) > e p(x) for all x G M 3 , i/ e = .B^ 1 - E > 0. (1.12) 



iv 



' /2 p)(0) = ^(M0) + ^ 2 p(0)). (1.13) 



dr 2 / 3 



Remark 1.12. The results in (i) generalize to the case of molecules, 
where the continuity results for p and h hold in the complement of the 
set {R x , ... , R L } C 1R 3 (see Remark |Oj ). 

Remark 1.13. It is known that eigenfunctions obey (Kato's) Cusp 
Condition (see Kato ||), and similar properties hold for particle den- 
sities. For more recent results see Hoffmann-Ostenhof et al. [|J, ||, 
Hoffmann-Ostenhof and Seiler |3j. In the proof of Theorem |1.11| , (iv) 
we make use of the Cusp Condition for p, namely: 

pV ) = lim ^ - ^ = -Zp(0) and limpYr) = p70) (1.14) 

and also present a proof for it. 

Remark 1.14. Of course our results are only first steps in a thorough 
investigation of qualitative properties of the one-electron density. Here 
are some obvious open questions: 

(i) Is p(x) > for all x G M 3 ? We remark that this cannot be true 
in general, since it is false for some exited states of Hydrogen. 

(ii) Is p G C°°(R 3 \ {0}) or even C W (M 3 \ {0}) ? 

(iii) Is p smooth in [0, oo), in the sense that ( Jxp) (r) exists for r > 
for all kl 

(iv) Is 4-p(r) < for r > ? This is expected to be true for 
groundstate densities, but not known even for the bosonic case 
like Helium. Our results imply that 4-p{r) < for r < Rq 
for the bosonic case, where i?o depends on the constant C in 
Theorem |1.2| . Note that because of ( |1.9|) and (|1.12| ) we have 



Ap > for \x\ > Z/e, and so the Maximum Principle gives that 
j^p{r) < for r > Z/e. 



Remark 1.15. In the proof of Theorem |1.11| we obtain (see Proposi- 

9 a a \ 

dxi i ' dxi 2 ' 3ii,3 ' ' 



tion O): With Vi = («^-, ^-), the function 



ti(r) = / / \V-LTp(ru,x 2 , ■ ■ ■ ,%n)\ dx 2 ...dxNdw 

JS 2 JR3(iV-l) 

is continuous on [0, oo). 
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2. Proofs 

Throughout the proofs, we will denote by C generic constants. 

Crucial for our investigations is Corollary 8.36 in Gilbarg and Tru- 
dinger ||. We shall make use of this result several times and for con- 
venience we state it already here, adapted for our special case: 

Proposition 2.1. Let Q be a bounded domain in Mr 1 and suppose u G 
W 1,2 (Q) is a weak solution of Au + Y^=i bjDjU + Wu = g in Q, where 
bj,W,g G L°°(n). Then u G C 1 '"^) for all a G (0,1) and for any 
domain Q' , Q' C Q we have 

\u\ C h>*(Q>) < C( sup \u\ + sup \g\) 
n n 

for C = C(n, M, dist(0', dtt)), with 

max {1, \\bj\\ L oo {n) , \\W\\ L °°(n), IMU°°(n)} < M. 

j=l,... ,n 

Thereby 

,, ,, „ \Vu(x) - Vu(y)\ 

\u\c^{n') = ||w||i,«(n') + || vu||i,oo(n/) + sup - . 

x,y£Q',x^y \% V\ 



Proof of Theorem |1.2| and Proposition |L5 . 



Let the function F be as in (|1.6|) and define the function F% by 

N g i 

Fi(xi,... ,x N ) = J2 ~ -^\j\ x j\ 2 + l + 1\/\ x 3 ~ x k\ 2 + 1 - 

3=1 l<j<k<N 

(2.1) 

A computation shows that 

AF = V, (2.2) 

\\F - -Fi||l°°(r3jv) , ||V(F - Fi)|| L oo( M 3jv) , II^FiHi 



<C(J3,N,Z), |/3|>1. (2.3) 

Make the 'Ansatz' ijj = e F ~ Fl ^i. Using (H - E)ip = and Q we get 
that ip\ satisfies the equation 

AV»i + 2V(F - Fx) ■ V^i + (|V(F - F x )\ 2 - AF 1 + E)ip x = 0. (2.4) 



Due to ( |2.3|) the coefficients in (|2.4j) are bounded in Mr . Then Propo- 
sition |2TT] implies that ip\ is C 1,a for all a G (0, 1), in any ball -B(x, R) C 
R 3N , and 

IV , i|c 1 ."(B(x,fl)) < C sup \4>i{y)\ (2.5) 

yeB(x,2fl) 

with C depending on R but not on x. Since 

|VV(y)| < |V(F - F 1 )\ |V(y)| + |e F - Fl V^i(y)| 
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we obtain, via (12.31) and (12.51), 



sup |W>(y)| < C{ sup |V(y)|+ sup |VVi(y)|) 

yeB(x,R) ySS(x,_R) yS-B(x,_R) 

<C(sup |V(y)|+ sup |Vi(y)|) < C sup |^(y) 

yeB(x,2i?) ye_B(x,2i?) yeB(x,2R) 



with C = C(R). This proves Theorem |L2 . 

Proposition p~~5| follows from ip = e F ~ Fl ipi, tpi G C^(B(x., R)), since 
e~ Fl is smooth. □ 



Proof of Theorem |1 . 11 . 

Multiplying the equation (H — E)ip = with ip and integrating over 
x 2 , ■ ■ ■ ,xn gives the equation 

f ipA\ip dx 2 ■ ■ ■ dxN + -j — r p{x\) = 

N r z 
= / / i^i — A,- — i — dx2 ■ ■ ■ dxN 

+ j -, - ip 2 dx 2 ■ ■ ■ dx^ 

= [ ^(H"- 1 - E)ipdx 2 --- dx N 

- r i 

+ >/ i r^p 2 dx 2 - ■■ dx N (2.6) 

where H N ~ l is the (iV — 1) - electron operator defined in ( |1.8|) . Since 
A^ 2 ) = 2|Vi^| 2 + 2^A!^ and / Ai(^ 2 )(xi, x') dx 1 = A 1 p in the 
distributional sense, we get that 

^Aip(xi) = / \Viip\ 2 dx 2 ■ ■ ■ dx N 

+ / i)A 1 ipdx 2 ■ ■ ■ dx N (2.7) 
which, together with (|2.(j|), gives the equation (r = |x|) 

J Ap(x) + -p( x )= ( V t^" -1 — E)ipdx 2 ■ ■ ■ dx N 

+ 7 / -j tjj 2 dx 2 ■ ■ ■ dxN 

JZ£ Jr3(n-i) |Xi - Xj\ 

+ / |Vi^| 2 rfx 2 ■ • ■ dx N = h(x), (2.8) 

jR3(JV-i) 

hence we obtain ( |1.9| ). Integration of ( |1.9| ) over 8 2 yields ( |1.10| ). 
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The proof of the regularity properties of the functions h and h are 
rather technical, and therefore postponed to the next section. 

We now verify the regularity properties of the functions p and p 
under the assumption that the regularity properties of h and h stated in 



Theorem |1 . 1 1| have been shown. Define the function p by the equation 

(r = |x|) 

p(x) = e- Zr (p(0)+p(x)). 
Then p = e Zr p — p(0), p(0) = 0, and (|2.8|) implies that 



(2.9) 



x 



Ap - 2Z - ■ Vp + Z A p = 2he Zr - Z 2 p{0). 



(2.10) 



be and since P e W^, 
also p G Wfa*. Therefore Proposition 10 leads to p G C!' a , for all 



Since h G L™ c , all coefficients of ( j2.10| ) are L 



a e (0, 1). Due to Q2 
Now consider 

Ap 



,PGC 



0,1/ 



) follows. 



2Z 

— p + 2h = g 

r 



in 



\{0}- 



(2.11) 



Since /i G C Q (1R 3 \ {0}), for all a G (0, 1) and due to the above, p/r G 
C a {R 3 \ {0}), for all a G (0, 1), we have 

#GC a (M 3 \{0}) for allae (0,1). 



(2.12) 



From ( [2. 11| ) and ( 2.12 ) we obtain from regularity theory for the Poisson 
equation that p G C 2,a (M. 3 \ {0}) (see e. g. Gilbarg and Trudinger 
Thm. 4.3 and 4.6] or Lieb and Loss |10|, Thm. 10.3]). 

We proceed analogously for p: integrating ( 2 . 1 1| ) over S 2 , we get the 
equation 



Ap 



with 



2Z. 



h(r) 



2h 



in 



\{0} 



h(ruo) duj. 



(2.13) 



(2.14) 



Since the R.H.S. of fl213D is in C Q (R 3 \ {0}), we obtain that p as a 
(radially symmetric) function in M 3 is C 2,a away from the origin, and 
therefore p : R + -> E satisfies p G C 2 ' a ((0, oo)). 

That p G C 2 ([0, oo)) is shown in the proof of (iv). 

Next we prove (iii). To prove the bound flnTJ), let E^' 1 be the 
groundstate energy for the operator H N ~ l . From Remark [TTS] and the 
Variational Principle we get that for almost all x% G R 3 , 

/ ip(H N - 1 -E)ipdx 2 --- dx N > (Eq' 1 — E)p(xi), 



and so h(x) > e p(x) with e = E { 



N-l 

o 



E > 0. 
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As for the bound note that due to the operator inequality 

—A—j3/r > —f3 2 /4 (true in dimension 3) and the translation invariance 
of — A we have, for almost all x^ G M 3 (fixed), k G {1, . . . , N}, k ^ j, 

I :\ip\ 2 dxj< j \V j ip\ 2 dx j + - j \i!\ 2 dxj. 

JR3 \Xj — Xk\ JR3 4 J R 3 

In this way, using (|2.7|) , 

h(x)<C( | Vtp\ 2 dx2 ■ ■ ■ dxN + / \tp\ 2 dx2 ■ ■ ■ dxN) ■ 

(2.15) 



Due to Theorem [L2] and a subsolution estimate (see Simon [[[J, The- 
orem C.1.2.]) we get, with x = • • • ,Xn) = (xi,x r ), x' G R 3 ^ -1 ) 
and xn the characteristic function of the set Q: 



|W(x)| 2 <C sup my)\ 2 <c[ my)\ 2 dy 

y€B(x,R) JyeB(x.,2R) 

= C XB(x,2/?)(y)|^(y)| 2 o?y- 

Using this, and that for fixed x, y G M. 3N : 



we get, by Fubini, 



1 if|x-y|<2 J R, 
otherwise 



/ \V^(x u x')\ 2 dx' 

<c f (f XB (x ,2R)(yM(y)\ 2 d y) & 

= C f my)\ 2 ( [ XB(y,2R)(x) dx') dy. (2.16) 

JR3N \ Jm.Hn-1) / 

Note that with z = x — y we have 

XB(y,2fl)(x) = XB{y,2R)(z + y) = XB(0,2H) (z) 

and so (with r' = |z'| and u/ = z'/r') 

XB(y,2R)((x 1 ,x'))dx' = / Xi?(0,2R) ( («1 , z') ) 

(AT-1) JR3(JV-1) 



dz' = Xb(o,2R)(zi) / / r' 3(7V ~ 1)_1 dr' gL/ 

(*!,*') [<2fl ,/s3(iV-l)-l JO 

= C7(iV)(4i? 2 -|^ 1 | 2 ) 3 ^ 1 )/ 2 XB{0i2 , ? , ) (^ 1 ) 

< C(N)R^ N -^ XB{xim ( yi ). (2.17) 
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From Q2.16Q and (|2.17|) we get 

\V^(x u x')\ 2 dx' 



JV-1) 

< C(R) 



\x!~ yi \<2R 



j 



m yi ,y')\'dy'd yi = C(R) / p( yi )dy 



B(x!,2R) 



(2.18) 



Combining (gJ^) and ( g3g ) proves (|TTTT] ). 

We now prove (iv). We first prove Kato's Cusp Condition ( |1.14| ) for 
the function p: 



7(0) = lim 

r|0 r 



-Zp(0) and limp'(r) = p'(O). 

r|0 



First, define the function p by the equation (see also (|2.9| )) 



p(r) 



(p(Q)+p(r)). 



(2.19) 



Note that /i(0) = 0. Then, using ( |1.10| ), p satisfies the equation 



x 



Ap - 2Z - ■ Vp + Z 2 p = 2he Zr - Z 2 p(0), 
r 



and, since h is continuous, Proposition |2.1| gives that (radially 
symmetric) function in M 3 , is C 1,a in a neighbourhood of the origin. In 
particular, lim r | p'(r) =p'(0). Since (see (|2.19|) ) 

p» = -Zp(0) + e- Zr p'{r) 



this means that 

lim p'(r) = -Zp(0) + limp'(r) = p'(0) - Zp(0). 

rlO rlO 

From ( |2.19| ) we also get that 

gOj = e zr P(r) - g0) + e Zr - 1 

This, together with ( [2.2 1| ) and p(0) = 0, implies that 

p[r) 



(2.20) 
(2.21) 



lim p'(r) = p'(0) - Zp(0) = lim !-±-L - p(Q) li m 

r|Q r|0 r r|0 r 

rlO r 

and by (^20|) and ( PT22|) , p'(0) = -Zp(0). This proves ( gig) . 
Next, define the function r/ by the equations 



p{r) 



-Zr 



(p(0)+(3r 2 + v(r)), 



(2.22) 



(2.23) 



1 ,~. . z- 
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Then 77(0) = 0, and due to ( PH) , 5f(0) = 0. Together with this 
gives 

a~/\ fid & \ ~ 
A V{r) = — jrr tt )v 
V r z or or J 

~ Z 2 

= 2 [e Zr h(r) - —p{r)e Zr - 3/3 + 2Z(3r + Zrf(r)] 

= 2G{r). (2.24) 

From the foregoing regularity properties of p, in particular ( P-.14|) and 
fl2.23|) , we obtain that rf G C°([0,oo)). From this, together with the 
regularity properties of h shown in Section |3|, we conclude that G E 
C°([0,oo)) and 

~ Z 2 

G{r) -> h(0) - ~ 3(3 = as r ^ °' ( 2 ' 25 ) 

From flESEp and (|2~24l) we get that 

7y( r ) =2 / \ [ G(s)s 2 dsdt, 
Jo t Jo 



and 

rj'(r) - rf(0) 2 



r 

so that 



hm! A I )- z A0) = 2 

r|o r 3 



Then by ( ggg) p"(0) exists, and 



r^) = z 2 m+w = \{h^) + z 2 m)- 



This verifies ( p..!3[) . 

Furthermore, by ( |2.24j ), 



and so 



2 

ArjCr) = T]"(r) + -r/(r) = 2G(r) 
r 



r7'(r) = 2G(r) - -rf{r) = 2G(r) - 2 ^ 
r 



since r/(0) = 0. This implies, by ( |2.26|) , 

hW'(r) = 2(limG(r) - lim ( ?±LllW ) ) = ff'(Q) = 0, 

r|0 W V r|0 r|0 \ r / y 

so that due to (|2.23|) p"(r) is continuous at r = 0. Hence formula 
( |1.13 ) follows from ( |2.23| ) and rj" (0) = 0. This finishes the proof of 



Theorem II. Ill □ 
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3. Regularity of h and h 



In this section we prove the statements in Theorem |1 . 1 1| on the regu- 
larity of the functions h and h. More precisely, we prove the following: 



Proposition 3.1. Let ip satisfy (|1.3| ) and let h be as in (|2T 
h(x) = / | Vi/j\ 2 dx2 ■ ■ ■ dxN 

7]R3(JV-1) 

N 



E 

i=2 



R 3(JV-1) \Xj\ 



?/> 2 dx2 ■ ■ ■ dXN 



+ E 

l<j<k<N 



(jv-i) — X fc 



?/> 2 da;2 • • ■ dxN, 



(3.1) 



and /i as in Q2.14|) : 



h(r) = / h(ruj) dev. 



Then h G C£ C (M 3 \ {0}) and h G C°([0,oo)) n C£ c ((0, oo)) for all 
a G (0,1). 

Remark 3.2. From the proof of Proposition |3.1| follows the continuity 
of the function t\ in Remark 1.15 . 

Proof. For convenience, we shall often write J = J^n-i) ■ Let 
Ji(x) = J | Vip\ 2 dx2 ■ ■ ■ dxN, 

N 

•m*) = E/ 



rip dx 2 ■ ■ ■ dx 



3=2 



E 



-ip 2 dx2 ■ ■ ■ dx 



N ■ 



(3.2) 



For the proof of the regularity of the functions J\, J2, J3 we shall make 
use of the following lemmas. The proof of the first lemma is trivial, 
using \(x h x 2 , • • • ,x N )\ > \(x 2 , • • • ,x N )\. 

Lemma 3.3. Let a G (0, 1) and xq el 3 . Assume that the real func- 
tion G = G(xi, . . . ,xn) satisfies: For all R > there exists constants 
C, 7 such that 

\G(x, x 2 ,..., x N ) - G(y, x 2 ,..., x N ) 



sup 

x,y£B(xo,R) 



< 



\x - y\ 



C exp(-7|(x ,x 2 ,... ,x N )\) for all (x , x 2 , . . . , x N ) G 



3^ 



(3.3) 
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Then the function 

7](x) = / G(x, x 2 , ■ ■ ■ , xn) dx 2 ■ ■ ■ oIxn 

JR3(JV-1) 

is in C? 0C (R 3 ). 

We next prove the following lemma. 

Lemma 3.4. Let a G (0,1). Assume that the real valued function 
K = K(x\, . . . , xjv) satisfies (|3.3|) and that there exists constants C, 7 
such that 

\K(x 1 , . . .,x N )\ 

<Cexp(-7|(xi,... ,2^)1) for all (x 1} . . . ,xjv) eR 3N . (3.4) 

For all j, k £ {1, ... , iV}, j 7^ fc, i/je function 

C(xx) = j 1 rif^i, • • • ,a;jv) f&r 2 - • -dx N 

y K 3(iv-i) |a?j — Xfc| 

zs zn C£ C (M 3 ). 

For j > 2 i/ie function 

/i(xi) = / -j — r K(xi, . . . ,xn) dx 2 - • ■ dxN 



is in C£ C (R 3 ) 



Proof. Assume first that j ^ 1 ^ k. Let x,y E B(xq, R), then by 
|C(^) — C(2/)l 



< 

x — y\ a 

1 \K(x,x 2 , . . . ,x N ) - K(y,x 2 , . . . ,x N ) 
\xj - x k \ \x-y\ a 



dx 2 ■ ■ ■ dx 



N 



< C -. r exp ( - 'y\(x ,x 2 , . . . ,x N )\) dx 2 ■ ■ -dx N . 



By equivalence of norms in R 3N there is a constant Co such that 

N . 

c( i J / exp(-7Co|^|) dx^j x 



IC(*)-C(y)l < 



F - y 



/• 1 

x / -j j- exp(— 7Co(|xj| + \xk\)) dxj dxu < C, x,y E B(xq, R). 

The last inequality is an application of the following inequality (with 
n = 3, \ = l,p = r = 6/5): (see Lieb and Loss [TO, Theorem 4.3]) 



14 
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Hardy-Littlewood-Sobolev Inequality: Letp,r > 1 andO < X < n 
with 1/p + X/n + l/r = 2. Let f e L p {R n ) and h G L r (R n ). Then there 
exists a sharp constant C(n,X,p), independent of f and h, such that 



f(x) \x — y\ x h(y) dx dy 



n .nan 



< C(n,X,p)\\f\\ p \\h\\r- 



This proves Lemma [3]4| (a) when j ^ 1 ^ k. 

Assume now that j = 1. We assume without loss that k = 2. Then, 
with xi,Xi G B(x , R), 

|C(*i)-C(5i) 



< 



Xi — Xi\ 

1 



+ 



\K(x 1 ,x 2 , ... , x N ) - K(x u x 2 , ■■■ ,x N ) 

Xi — x 2 \ 

1 1 



dx 2 ■ ■ ■ dxN 



X\ — x 2 \ 
(A) + (B). 



\x 1 - x 2 \ 



\x\ — x~i\~ a \K(xi, x 2 , . . . , xn)\ dx 2 ■ ■ ■ dXjq 



For (A), using (|3.3|) and equivalence of norms in 



we get 



(A) < C ( rexp ( - 7c (|x | + \x 2 \ H h \xn\)) dx 2 ■ ■ ■ dx 

J \xi-x 2 \ 

exp(-7c |x 2 |) dx 2 < C(x 0: R). 



< C 



\X\ - x 2 \ 



As for (B), we apply the following inequality; for the convenience of 
the reader, we give the proof (borrowed from Lieb and Loss [0, (3) p. 
225]). 

For a G (0,1): 

1 1 



x — z 
x — z\~ 



1-a 



\y - A 



+ \y-z 



x-y\- a < 
x ~ a for all x, y, z G M 3 . 



(3.5) 



Proof of (|3.5|). By Holder's inequality we have, for b > 1, a G (0, 1) 



1-b- 1 



r 2 dt < 



dt 



t- 2 /^ dt 



l-a 



<(b-l] 



Substituting b/a for b, with a > 0, this gives 



I6" 1 



.-i I 



< \b — a\ a max{a 



— 1— a v—l—a 



}• 



So, for x,y,z G M 3 , using max{s, t} < s + 1 and the triangle inequality 
in IR 3 , we have 



\x — z\ 



\y - z 



< \x-y\ a { 



x — z 



-l-a 



+ \y- Z\ 



-l-a 



}■ 



□ 
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In this way, by Q3.4|) and equivalence of norms in M. 3N , 

I f \K{x u x*. x N )\ + \K(x 1 ,x 2 ,. x N )\ \ 
J V x 1 -x 2 1+a \x l -x 2 \ 1+a J 



N f 

I / ex P( — 7 c okj|) dxj 

7=1 ^ 



x y R3 ( - — ztttt^ + - — znrr^ ) exp(- 7 co|x 2 |) dx 2 
<C(x ,R), 



R3 -x 2 \ 1+Q \x 1 -x 2 \ 1+a 



since xi, x\ G B(xq, R). This finishes the proof of Lemma |3.4| (a). 
The proof of (b) is similar to that of (a) so we omit the details. □ 

The proof of the following fact is straightforward: 

There exist constants C = C( 7 , R) and 7 = 7(7) such that 

exp(-7|(x, ... ,x N )\) <C exp(- 7 |(x , ... ,x N )\) (3.6) 

for all x G B(xq, R). 



Using this and Lemma |3.4| , we shall prove the following lemma 
on the regularity of the functions J\, J 2 and J3 from ( |3.2| ). 



Lemma 3.5. Let J 1; J 2 and J3 be as in (|3.2| ). Then 

(i) J 2 ,J 3 eC? oc (R 3 ) for alia e (0,1). 

(ii) Ji G ^(R 3 \ {0}) /or a// a G (0, 1). 

Herefrom follow the regularity properties of the function h stated in 



Proposition |3TT 



Proof of Lemma |3.5| (i) . 

Firstly, by Theorem |1.2j and Remark |L8 



|V(x)| , |V^(x)| < Cexp(- 7 |x|) , x G R , (3.7) 

which gives (|3.4|) for K = ip 2 . 

Next we verify that for G = ip 2 , (|3.3| ) is fulfilled. Then Lemma 3^ 
can be applied with K = tp 2 , and the Holder-continuity of J 2 and J3 
follows. 

Given xo G R 3 ,i? > 0,a G (0,1), and x, y G B(xq,R). Using that 
(see e. g. Maly and Ziemer JLT], Theorem 1.41]) (here, (x, x 2 , . . . , xjy) = 
(x,x'), x' G R 3 ^- 1 )) 

V> 2 (x, x') - V 2 (y, x') = y 7^ [V> 2 (sx + (1 - s)y, x')] ds 



[ [V 1 (ij 2 )(sx + (l-s)y,x')]-(x-y)ds 
Jo 



(3.8) 
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and that sx + (1 — s)y G B(x , R) for all s G [0, 1] we get, with (|3 
and (KB 



\ip 2 (x,x') -ip 2 (y } x')\ 
\x — y\ a 

<2\x-y\ 1 ~ a [ \Viif)(sx + {l- s)y,x')\-\if)(sx + (l- s)y,x')\ds 



o 



<2(2R) 1 ~ a / C exp(-2 7 |(sx + (1 - s)y,x')\)ds 
Jo 

< C exp(-7|(x , • • • ,x N )\), (3.9) 

so (|3~3| ) follows for a G (0, 1). 
This proves (i) of Lemma |3.5| . 

To prove (ii), we write if) as in the proof of Theorem |1.2| : if> = 
e F - Fl i() 1 , with F and F 1 as in Q) and (g^). Then 

, lW = j|V^< = /|VF|V^ + /| W ^ 
-2 y (VF- VF 1 )^ 2 dx'+ / e^-^lWil 2 ^' 

+ 2 y (VF- V^i)e 2(F - Fl Vi^'-2 y (VFx-V^Oe^-^Vi ^' 
= + h{x) + 7 3 (x) + I 4 (x) + I 5 (x) + I 6 (x). (3.10) 
Using the idea from (|3.8|) twice (on |V.Fi| 2 and ?/> 2 , respectively), the 



estimates (|3"77|), Q, and (PQ, we have, with x G M 3 , i? > 

0, a G (0, 1), and x, y G B(xq, R): 



SJF 1 \ 2 if) 2 (x,x') - \VF l \ 2 if) 2 {y,x')\ 
\x — y\ a 
\\\7F 1 (x,x>)\ 2 -\VF 1 (y,x')\ 2 \ 



< 



t\\2 



,„ 1^(^)1 

F — 2/1 

,x ,o l'0 2 (^ ) — ?/> 2 (y, x'x)! 

+ivf 1 ( 1/ ,xoi 2 1 ^ Jl 

< C |a; - y| 1 - Q ||V(|VF 1 | 2 )|| oo exp(-2 7 |(x,a; 2 , . . . ,0^)1) 

+ 2C\x - 2 /| 1 - Q |||VF 1 | 2 || oo exp(-7|(x , ...,x N )\) 

< C exp(-7(x , • • • ,x N )\). 



By Lemma |373| , with G = |V-Fi| if) (xx, ■ ■ ■ , xjy), this implies that J 2 G 

cg c (K 3 ). 

Using the same ingredients, writing V^i = V(e Fl F ip), gives (|3 
and (U) with 
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and 



G = K = (VFx ■ V^i)e 2(F - Fl Vi, 



and so by Lemma I4, J 6 e (^(IR 3 ). 

The remaining terms are those involving the function F, namely I±, 
h and J 5 . 

Note that 



VF(xi, . . . ,x N ) = - — ( J 1 



+ 4 



and so 



, , N N N-l 

1 / Xi - Xfc ^ Xj - x k 

A \ < ' Iti T*7. I ' ' f ^ It - TiJ ' ' ' J 



x N - x k 



k=2 \ Xl -X k \ k=l ^. \ Xj ~X k \ k=i \X N - X k \ 



(3.11) 



\VF{x!,... ,x N )f 



NZ 2 Z Xj Xj — x k 

~~ 4 ~8 ^ 

N 



+ 



1 E 



X j | | X j X | 

X j X X n X i 



16 •,, TTI-w \ X 3 ~ X k\ \ X j ~ X l\ 



In this way, 

nz 2 r 

h(x\) = — - — / ip 2 dx2 ■ ■ ■ dx 



N 

N r 

-t E / rV i^'.V ^-^ 

1 JV _ _ 

1 \ ^ Xj X k Xj X[ 2 7 , 

+ — > 1 • r V> ^2 • • • dxjy 

z l * 

POi) - ^ E + 7^ E ^.wfri)- ( 3 - 12 ) 



NZ Z , , Z 



16 

j,k=l,k^j j,k,l=l,k^j,l^j 



Note that i/ 3 - (fc)J = u jt i ik . 

Using the ideas above, Lemma |3.3| implies that the following func- 
tions from ( |3.12| ) (with the mentioned choices of G satisfying ( |3.3| )) are 
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all inC? oc (M 3 ): 



p ■ 


G = 












G = 


Xj 

\ x j 


Xj 

\ x j 


- x k 

- x k 




Vj,k,k> 3 7^ h : 


G = 


Xj 

Fj 


- X k 

- x k \ 


Xj 

\ x j 


- x k 

- x k 




G = 


Xj 


- x k 


Xj 


- Xi 


\ x j 


- x k \ 


\ x j 


- Xi\ 



Vjhu j,k,l^l,l^j^k: G = - 1 -i) 



Likewise, Lemma |3.4j implies (with the mentioned choices of G = K 
satisfying ( |3.3| ) and (|3.4|) ) that the following functions from ( |3.12| ) are 
all inCT oc (M 3 ): 



K jA , jvi: g = k= j y, ^ 

\x 



Xj ' yJJC j X 'Y J .) 

'il 



j,l =£1,3^1: G = K = —t — | 1 V, 

From the decomposition of l\ in ( |3 . 1 2|) we are left with 

«i, fc (si)= / A" , Xl ~ Xk , ^ 2 dx' , fc = 2,...,JV, (3.13) 
7 |x x | \xi-x k \ 

and 

"i,fc,t(si) = [ l Xl ~ X \ ■ Xl ~ Xl tfdx* , A;,/G{2,...,iV},M^ 

(3.14) 



Note that 



X\ X\ X k 2 i i 

■ dx 2 • • • ax at 



Fi 



Xl\ J \xi-x k 



■ (xi • (xi — x k ) tp 2 ) dx 2 ■ ■ ■ c?x at. 



The function l/|xi| is smooth for x\ ^ and therefore in C 1 " c (M 3 \{0}). 
The function x\ ■ (x\—x k ) ip 2 satisfies ( |3.3|) and (|3.4j) (by the same ideas 
as above), so Lemma |3.4| (a) implies that the function 



Fi _ x k 



■ (xi • (xi — Xk) 4> 2 ) dx 2 • • • dx 



N 



is in C£ C (M 3 ). The functions in ( g7LjD are therefore in C? oc {R 3 \ {0}). 

As for the functions in ( 3.14Q , these are all in Cj" c (R ), which can 
be seen by applying the previous ideas, in particular Lemma |3~5| , 
O and (0). 

This proves that h G C£ C (M 3 \ {0}). 
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As for J 3 (see (ETTOp and (ETT1I)), with V = (Vi, . . . , Vjv), 

N „ 

7 3 (x) = Z^ / (p- - V j F 1 Wdx 2 ---dx N 



3=1 



-o £ /(i^l-V^y^..-^. (3.15) 

3,k=l,3i=k 



The terms in the first sum with j ^ 1 are in Ci" c (IR 3 ), due to Lemma 3^ 
(b), with G = K = (xj-VjFt) ^ satisfying (pD and (B- ( To see this, 
use the previous ideas; to apply the idea from (|3.8|) to Vj.Fi we use that 
Fi is smooth). The terms in the second sum in ( |3.15| ) are all in Ci" c (IR 3 ), 
due to Lemma |3]4| (a), applied with G = K = ((xj — Xk) ■ Vj.Fi) ip 2 - 
The term with j = 1 is in Ci" c (R 3 \ {0}). This can be seen by following 
the ideas in the proof of the regularity properties of the function in 
(|3.13| ), now using Lemma |3]3| with G — (xi • ViFi)^ 2 . 



The statements and proofs are similar for 

N „ 

h{x) = / (A • V^i)e 2(F " Fl Vi dx 2 --- dx 

._, J ^ \ x j ' 



3=1 

N 



1 y, I (^-^.V^ 1 )e^ F - F ^ l dx 2 ---dx N . 

(3.16) 



2 

j,k=l,jj^k 



That is, the functions in the first sum in ( |3.16j ) with j > 2 and those 



in the second sum are all in C£; C (R ), whereas the function in the first 
sum with j = 1 is only in Ci" c (IR 3 \ {0}). To prove this we use the 
inequality (with x = (x±, . . . , £jv)): 

\iPi\c^(b(x,r/2)) < C sup \if>i(y)\ < C exp(-7|(x 1 , . . . ,x N )\). 

ye(B(x,R)) 



This inequality follows from (|2.3|) , (|2.5| ) and ( |3.7|) (remember that ^i = 
e Fl ~ F i)). 

This proves that I 5 G Ci" c (IR 3 \ {0}), and so finishes the proof that 
Ji G C£ C (M 3 \ {0}). (See (pT0|) ). This proves (ii) and therefore 
Lemma |3.5| . □ 

That /i G Ci" c ((0, oo)) is a consequence of the foregoing and of the 
following proposition: 

Proposition 3.6. Assume f G C£ C (M 3 \ {0}), a G (0, 1). Then f G 
Ci Q oc ((0,oo)) ; where 

7( r ) = / f( ruj ) duj - 

Js 2 
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Proof. Let r G (0, oo). For all x G A = {x G IR 3 | |x| = r}, choose 
i? = R(xq) and C = C(xo) such that 

sup \JM^1M < C(xo) . (3.17) 

x,seB(. *o)) F - S/r 

This is possible, since / G C^M 3 \ {0}). Then 

Ac |J B(x ,R(x )). 

Using compactness of A, choose x±, . . . ,x m € A such that 

A c \jB(x j ,R(x j )). 

Choose e G (0, r) such that 

m 

{y E M 3 \ r - e < \y\ < r + e} C [J B(xj, R(xj)). 

j=i 

Then, for all s,t G (r — e, r + e) and all u G § 2 there exists j G 
{1, . . . , m} such that sou, tu G -B(£j, R(xj)) and therefore by ( |3.17| ), 



|s — t| a I so; — tu; 

So with C = max{C(xi), . . . , C(x m )}, 

\f(su)-f(tLu)\ 



< C(x s ). 



< C, for all s, t G (r — e, r + e) and all a; G S . 
This implies that 

\7(s)-7(t)\ _ iy/H-/(M)^i 

This proves that / G C£ c ((0, oo)). □ 

To prove that h G C°([0, oo)), we apply the following: 
Proposition 3.7. Assume f G Cg c (M 3 ). Then J e C°([0,oo)), where 

7( r ) = / f( ruj ) duJ . 



Proof. The function / is continuous in R 3 , since it is in C^IR 3 ). Let 
r G [0, oo). Then 

lim/(su;) = f(ru) for all uj G S 2 . 
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Using the supremum of / on a sufficiently large compact set in M 3 as a 
dominant, Lebesque's Dominated Convergence Theorem gives us that 



lim 



f(su) duj = I f(ru) dw. 



□ 



Therefore / G C°([0,oo)). 

Recall the proof of the fact that h G C£ C (M 3 \ {0}). In fact, the only 
terms in the decomposition of h (see ( |3.1| ), ( |3~2|) , (|3.10| ), and (|3.11| )) 
that are only in C^IR 3 \ {0}) and not in C^M 3 ) are the functions 



X\ - x k 



X\\ \X\ 
X\ 

Nil 



Xk\ 



i/j 2 dx2 ■ ■ ■ dxN , k = 2,...,N, 



ViFijip 2 dx2 ■ ■ ■ dxN, 
I (^-Vi^e^Vi^ 



dxN- 



(3.18) 



Comparing ( |3.10| ), (|3.13|) , (|3.15|) and (|3.16|) , it can be seen that all the 
terms in (|3.18|) stem from the function Ji, namely from Ji, J3, and I5. 

All other terms in the decomposition of h are in C^R 3 ). When inte- 
grating them over § 2 , we get something continuous in [0, 00), according 
to Proposition [T7| above. 

For the terms in ( 3.18|) we note that they are all of the form 

[ P- ■ K( Xl ,x') dx' = P- ■ f K(x u x') dx'. (3.19) 
J Fi| \xi\ J 

In each case, we have 

~L(xi) = (Li(xi), L 2 (xi), L 3 (xi)) = J ~K(xi,x')dx' ,Lj e Cj 

(3.20) 

To see this, apply Lemma to each of the coordinate functions Lj, 
j = 1, 2, 3. The integrands are easily seen to satisfy ( p.3| ) in each case, 
by the previous ideas. To get continuity in [0, 00) of the functions in 
( |3.18| ) we use ( |3.19| ) and ( |3.20D , and the following lemma: 

Proposition 3.8. Assume f = (fi, f 2 , fs), fj e C^ QC (R 3 ). Then f G 
C°([0,cx))), where 



r oc (K 3 ) 



fir) 



[uj ■ f(rtu)) duj. 



Proof. The same as for Proposition |3.7| , noting that for all r G [0, 00) 
and fixed uj G S 2 : 

limu; • f(su) = uj ■ f(ruj). 



□ 
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This holds even in the case r = 0, for which 

lim / u ■ f(suj) du = I uo ■ f(0) duj = 0. 

This proves that the functions in (|3.18| ) are in C°([0, oo)). Therefore 
h G C°([0, oo)), which finishes the proof of Proposition |Q| . □ 

Acknowledgement: The authors wish to thank G. Friesecke and 
S. Fournais for stimulating discussions. 
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